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Chemically non-equilibrated quark-antiquark matter is studied within 
the Nambu- Jona-Lasinio model. The equations of state of non-strange 



' (q = u, d) and strange (q = s) qq systems are calculated in the mean-field 

ON ' 

approximation. The existence of metastable bound states with zero pressure 
is predicted at finite densities and temperatures T< 50 MeV. It is shown 

00 " 

ON ' that the minimum energy per particle occurs for symmetric systems, with 

Q_i' equal densities of quarks and antiquarks. At T = these metastable states 

have quark number densities of about 0.5 fm~ 3 for q = u, d and of 1 fm~ 3 for 



q = s. A first order chiral phase transition is found at finite densities and tem- 



^ ' peratures. The critical temperature for this phase transition is approximately 



75 MeV (90 MeV) for the non-strange (strange) baryon-free quark-antiquark 
matter. For realistic choices of parameters, the model does not predict a phase 
transition in chemically equilibrated systems. Possible decay channels of the 
metastable qq droplets and their signatures in relativistic heavy-ion collisions 
are discussed. 

PACS numbers: 12.38.Mh, 11.30.Qc, 12.39.Fe, 25.75.-q 



1 



I. INTRODUCTION 



Studying the equation of state (EOS) of strongly interacting matter is necessary for 
understanding the evolution of the early universe, properties of neutron stars and dynamics 
of heavy-ion collisions. First of all, one should know the degrees of freedom which are most 
relevant at a given energy density. When studying the EOS it is also important to identify 
the regions of possible phase transitions and to find all stable and metastable states which 
correspond to extrema of a thermodynamic potential. Presumably, lattice QCD calculations 
can be used to study these questions from "first principles". At present, however, this 
approach gives reliable results only for a baryon-free matter in thermodynamic equilibrium. 

Therefore, one should use effective theories when dealing with a wider class of multipar- 
ticle systems. The Nambu — Jona-Lasinio (NJL) model is one of such effective theo- 
ries which proved to be rather successful in describing the ground states of light hadrons 
in vacuum. The chiral symmetry of the strong interactions, most important for the low 
energy hadron physics, is explicitly implemented in this model. Although the model is non- 
renormalizable and does not contain gluons [], it seems to be a reasonable approximation to 
QCD at scales intermediate between asymptotic freedom and confinement. 

In this paper the NJL model is applied to investigate the properties of chemically non- 
equilibrated quark-antiquark (qq) matter. The study of possible bound states and phase 
transitions is one of the main goals of the present work. The search for new states of strongly 
interacting systems, essentially different from normal nuclear matter, has already quite a long 
history. Using the linear a-model, Lee and Wick predicted |§ a metastable superdense state 
of nucleonic matter at T = 0. This problem was further studied in Ref. |4j] by including the 
repulsive vector interaction. The liquid-gas as well as the chiral phase transitions of hadronic 



Lattice calculations show that at temperatures T<150MeV the gluons acquire a large effective 
mass due to the color screening effects. In the present work, dealing mostly with such moderate 
temperatures, we assume that gluonic degrees of freedom are suppressed. 
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matter have been investigated within the framework of the relativistic mean-field model in 
Refs. 0-0]. Since the chiral symmetry is not respected by this model, its predictions become 
questionable at high baryon densities and temperatures. Attempts to generalize the hadronic 
models by implementing the main symmetries of QCD have been made in Refs. p|-|T0[|. A 
chiral phase transition was found in nuclear matter only at very high baryon densities and 
temperatures when hadronic degrees of freedom are inappropriate. 

The existence of some exotic multiquark and multihadron states has been discussed by 
various authors. For example, "stranglets" (bound states of strange and light quarks) have 
been proposed in Ref. fLl| . The possibility of their formation in relativistic heavy-ion col- 
lisions has been considered in Ref. |fL2|| . Arguments that multipion systems may be bound 
have been given in Ref. [plj. The phase transition of quark matter into a color supercon- 
ducting state has been predicted recently fll^,|T5|| within the QCD-motivated models. The 



random matrix model has been used to study the phase diagram of quark matter in Ref. |T6 
Possible signatures of QCD phase transitions have been discussed in Ref. [[Hj . 

The possibility of multiquark bound states and chiral phase transitions was studied by 
many authors (see e.g. JT8|-pn|]) within the NJL model. Attempts to apply the same model 
to nucleonic matter were made in Refs. p2|-p^|] . Most of these studies have been done within 



the mean-field approximation. An attempt to go beyond this approximation, by including 
the mesonic degrees of freedom, has been made in Ref. |25fl . The results of these works, 
obtained with different sets of model parameters, often contradict to each other. A common 



shortcoming appearing in many papers is the omission pO|j2T[| or underestimation [|T8j of the 



repulsive vector interaction^. It was shown in Ref. Jl9| that this interaction becomes rather 



important at large baryon densities. According to Ref. [|19j, at realistic values of the vector 
coupling constant the NJL model does not predict a first order phase transitions in quark 
matter. 



The vector interaction was not considered in the original version of the model 
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All these results have been obtained by assuming that the conditions of statistical equilib- 
rium are perfectly fulfilled in the matter. This implies two kinds of equilibration. The first, 
kinetic or thermal equilibrium, means that the occupation numbers of quarks and antiquarks 
coincide with the equilibrium distribution functions characterized by some temperature and 
chemical potentials. The second, chemical equilibrium, requires that the certain relations 
between chemical potentials of different particles must be fulfilled. These relations deter- 
mine the equilibrium abundances of various species which can be reached at large times due 
to inelastic multiparticle interactions. First attempts to study non-equilibrated qq systems 
were made in Refs. |26|-p8| by using transport equations derived from the NJL model. 

Off-equilibrium effects should inevitably accompany the formation and evolution of quark 
matter in relativistic heavy-ion collisions. There are many theoretical as well as experimen- 
tal arguments in favour of large deviations from the chemical equilibrium, even in central 
collision of heaviest nuclei at cm. bombarding energies \/s > 10 GeV per nucleon (i.e. for 
SPS energies and higher). This is a consequence of a short time available for the interaction 
of primary nuclei and for subsequent equilibration of particles at such high energies. The 
degree of the chemical equilibration has been investigated p9| , |30| on the basis of the parton 
cascade model. Strong deviations from chemical equilibrium were predicted for light quarks 
at intermediate stages of a heavy-ion collision at RHIC and LHC energies. Large nonequi- 



librium effects at SPS energies were found within the UrQMD model pl |. As follows from 



the thermal model analysis of experimental data |H| , the hadronic matter is out of chemical 



equilibrium at late stages of relativistic nuclear collisions. The conclusion about a significant 



"overpopulation" of light quarks in 160 AGeV Pb+Pb collisions has been made in Ref. [|33 

Motivated by these findings, below we study the properties of chemically non- 
equilibrated qq matter which is however in the thermal equilibrium. Formally, we intro- 
duce two chemical potentials for quarks and antiquarks with a given flavour and assume 
that their values do not satisfy in general the conditions of chemical equilibrium. It will be 
shown below that the EOS of chemically non-equilibrated qq matter is non-trivial and its 
phase structure is much richer as compared with the equilibrated matter. 



In the next section we briefly formulate the SU(3) version of the NJL model used in the 
present paper. In Sec. 3 this model is applied to calculate the EOS of the non-strange and 
strange qq matter at arbitrary densities of quark and antiquarks. The numerical results are 
presented in Sec. 4. Sec. 5 is reserved for discussion and summary. 



II. FORMULATION OF THE MODEL 



We proceed from the SU(3) version of the NJL model suggested in Ref. Its inter- 
action part is given by the local coupling between the quark color currents. After the Fierz 
transformation the color singlet part of the Lagrangian may be written as: 



C = ijj (i d - mo) i> + G s 

3=0 
8 

-gvy: 

3=0 



% A 
2 



(1) 



where Ai, . . . , \% are the SU(3) Gell-Mann matrices in flavour space, Ao = J 2/3 7 and 
mo = diag(mo u , mod, mos) is the matrix of bare (current) quark masses. At mo = this 
Lagrangian is invariant with respect to the SU(3)l <8> SU(3)r chiral transformation. The 
term with mo leads to explicit breaking of chiral symmetry which is supposed to be small. 
The relation 



G,q — 2 G\ 



(2) 



between the scalar (Gs) and vector (Gy) coupling constants follows from the QCD motivated 
initial Lagrangian |34[ . 



In fact, the above model can be formulated in such a way that Gy/Gs is considered as 
an adjustable parameter, which may be chosen by fitting the observed masses of the vector 



mesons. Different calculations use values of Gy/Gs in the range 0.5-1 p4||35| . We would 
like to stress here that the choice Gy = 0, frequently made in the literature, is unrealistic for 
baryon-rich qq systems. In this case, the omission of the repulsive vector interaction results 
in an underestimation of pressure leading to incorrect results concerning the possibility of 



phase transitions and the existence of (meta)stable states (see also Ref. p3f). On the other 
hand, as will be seen below, thermodynamic functions of baryon-free systems are insensitive 
to G v . 

In the mean-field (Hartree) approximation, used in the present paper, only the scalar 
and vector terms of Eq. ([I]) survive and the Lagrangian is diagonal in the flavour space: 

£mfa = /] £>f , (3) 
f=u, d, s 

where 

Q Q 

C f = ^ f (iD - m f )i(> f - — < i>$t > 2 + ~y < ^ftv&t >2 • ( 4 ) 

Here ipf denote spinors for u, d, s quarks, = d^ + iGy < ip p^ipj >, angular brackets cor- 
respond to quantum-statistical averaging and the constituent quark mass m/ is determined 
by the gap equations 

rrif = m f - G s < ip f ipf > ■ (5) 

Due to the absence of a flavour-mixing interaction, all extensive thermodynamical functions 
(energy density, pressure etc.) are additive in quark flavour. Therefore, in this approximation 
light (/ = u, d) and strange (/ = s) quark systems can be studied separately. In the 
following we consider isospin-symmetric qq systems, assuming equal numbers of u and d 
quarks (antiquarks), and disregard the difference between m 0u and m 0c i . 

The model parameters mo/,G5,A (the ultraviolet cutoff in 3-dimensional momentum 



space) can be fixed by reproducing the empirical values |36| of quark condensates < uu >, 
< ss > and the observed values of 7r and K decay constants fx, fx Q in the vacuum. We 
choose the following values of the model parameters 

m 0u = mod = 7 MeV, m 0s = 132 MeV, G s = 24.5 GeV~ 2 , A = 0.59 GeV , (6) 



3 These decay constants are related to the bare quark masses and the condensate densities by the 
GOR relations |7|. 
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The calculation shows (see Sec. 3) that these parameters correspond to the values 

U = 93 MeV, < uu >= (-230 MeV) 3 , m™ c = 300 MeV , (7) 
f K = 90 MeV, < ss >= (-250 MeV) 3 , m™ c = 520 MeV , (8) 

where m™ c is the constituent quark mass in the vacuum. 

III. EQUATION OF STATE OF QUARK ANTIQUARK MATTER 

In the mean field approximation the Lagrangian is quadratic in the quark fields and the 



calculation of thermodynamic functions is straightforward. Below the method of Ref. [18 
is generalized for chemically non-equilibrated systems. The Coulomb and surface effects are 
disregarded and the characteristics of qq matter are assumed to be spatially homogeneous 
and time- independent. In this article we consider separately the systems with light (/ = u, d) 
and strange (/ = s) quarks. Unless stated otherwise, the flavour index / will be omitted. 
The main feature of a chemically non-equilibrated system is the existence of two chemical 
potentials, i.e. for quarks (p) and antiquarks (/Z) which are in general not strictly related 
to each other. If the equilibrium with respect to creation and annihilation of qq pairs were 
perfectly fulfilled, the above chemical potentials would satisfy the condition of chemical 
equilibrium, 

TJ=-/i. (9) 

At given temperature T and quark density p Eq. @ fixes the density of antiquarks p (see 
Fig. 8). Below we study a general case considering quark and antiquark densities and 
accordingly, their chemical potentials, as independent quantities. 

The generalized expression for the partition function Z of the thermally equilibrated 
qq system may be written as 

Z = Spe- 5 / T = e- Q / T , (10) 

where 



H = H - pN --pN . (11) 

Here H is the Hamiltonian and A" (N) is the quark (antiquark) number operator. Up to an 
arbitrary additive constant the thermodynamic potential Q is equal to — PV where P and V 
are respectively the pressure and the volume of the qq matter. Below the quantum-statistical 
averaging of any operator A is made in accordance with the expression: 

< A>=Sp(Ae~ 6 / T )/Z. (12) 

From the Lagrangian, Eq. (J|), one obtains the equation of motion for the quark field ip 

(i<y> t d lt -m-G v p v >f)ip = 0, (13) 

where m is the constituent quark mass and pv =< ip^ip > is the vector density. It can be 
expressed through the densities p =< N > /V and p =< N > /V as0: 

p v = p-p. (14) 



After the plane wave decomposition |39| of ip, ip , the quantum states of quarks and 



antiquarks can be specified by momenta p and discrete quantum numbers A (helicity and 
color). Let ap t \ (bp t \) and Op x (bp t \) be the destruction and creation operators of a quark 
(an antiquark) in the state p, A . Then the operators H, N and can be expressed through 
the linear combinations of Op A ap^\ and bp A bp t \ . For example 

77 = £ a p,x a P,A, ^ = b p,x (15) 

p,A V p,X 

In the quasiclassical approximation the sum over momenta may be replaced by the 
3-dimensional integral over momentum space: 

£ = 7^3 / rf V ( 16 ) 



p 



(27TY 



J The baryonic density ps of the qq system equals l/3py . 



8 



By using Eqs. (H)— © one can find the Hamiltonian density of the qq matter in the 
mean-field approximation. The resulting expression for the operator H reads as 

y = J2( E p- Vr) a p,A a f>,A + J2( E P - Vr) b p,x h P+ 



p,x 

~J2Ep + 

p,x 



p,x 

2 ry „2 



(m - mg) _ G y p y 
2 Go 



(17) 



where Ep = ^Jrri 1 + p 2 and //r, p R denote the "reduced" chemical potentials: 

{M R = jl- Gypy , 
JI R =JI+ Gypy . 



(18) 

(19) 



In Eq. (JTyD the gap equation 



m = m - G s ps 



(20) 



has been used. Here ps =< > is the scalar density of the qq system. It can be expressed 



as 



Ps =< IT ( a p.A a P,A + b p,\ bp,x-D> ■ 
p,x h P 



(21) 



Eqs. ([IT]) and (21) contain divergent terms. They originate from the negative energy levels 
of the Dirac sea [51,^3]. In the present paper these terms are regularized by introducing 
the three-dimensional momentum cutoff 9 (A — \p\) where 6 (x) = \ (1 + sgnx). The cutoff 
momentum A is considered as an adjustable model parameter. 

The quantum-statistical averaging of the operator ap~ x ap t \ gives the occupation number 
of quarks in the state p, A. Direct calculation shows that it coincides with the Fermi-Dirac 
distribution for an ideal gas with the chemical potential p R . 

Ep — hr 



< ap X ap^ >=n p 



cxp 



T 



+ 1 



(22) 



The analogous expression for the antiquark occupation number, np =< bp~ >x bp t \ >, is 
obtained by the replacement pn — > p R . At given densities of quarks and antiquarks, the 
chemical potentials p R and ~p R are determined from the normalization conditions 



P P 

where v is the spin-color degeneracy factor. For a qq system with single flavour 



(23) 



v = 2N C = 6 . 



(24) 



After introducing the ultraviolet cutoff, Eq. (PT|) takes the form 



fry r- 

Ps = v^2— np + n p -8(A-p) 
p h P 



(25) 



The explicit form of the gap equation is given by Eqs. (|20|), (p2|), (|25|). The physical vacuum 
(p = p = Pv = 0) corresponds to the limit rip = rip = 0. As discussed above, the pressure 
can be represented as 



T 

P = — In Z + P 
V 



(26) 



where the additive constant Pq is introduced to ensure that the vacuum pressure is zero, 
-Pvac = 0. Direct calculation gives the result 

GvPv 



P = Pk + 



B(m) . 



(27) 



Here the first ("kinetic") term in the r.h.s. coincides with the total pressure of ideal gases 
of quarks and antiquarks having constituent mass m and chemical potentials pr and ~p R : 



P k = vTJ2 In 
P I 



1 + exp 



Pr — Ep 
f 



+ In 



1 + exp 



P 



T 



v 



E TET ( n P + n P 



p 3E p 



(28) 



The third term in Eq. ( P7[ ) gives the "bag" part of the pressure. The explicit expression 
for B(m) can be written as 



B(m) = $(m vac ) - $(m) 



(29) 



Here m vac is the constituent quark mass in the vacuum and 

(m — mo) 2 



$( m ) = uJ2E p e(A-p) 



P 



2G f 



(30) 
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From Eqs. fl2"7|)— (p9|) one may see that P vac = —B (m vac ) = 0, i.e. the condition of vanishing 
vacuum pressure is automatically satisfied. The expression for $(m) can be calculated 
analytically, 

, . v A 4 f m\ (m — inn) 2 , . 

*( m » = 8^*(A)- L W^' < 31) 



where 



(x) = Ajdtt 2 Vt 2 + x 2 = (l + ^pj VTTx^ 
o ^ ' 



1 



1 + Vl + a; 2 



— ln^ ■ . (32) 



x 



As required by the thermodynamic consistency, the gap equation (^) should follow from 
the minimization of the thermodynamic potential Q with respect to m, 

fa d\ rriQ-m 

( d mP) T ,^ji = —q~ s Ps = 0. (33) 

Therefore, the quark mass is not an independent variable and the pressure can be considered 
as a function of T, fi and ~p only. The gap equation in vacuum is f\ 

$'(m vac )=0. (34) 

This equation was solved numerically with the values of model parameters from Eq. (^j). 
The calculations were performed separately for light (/ = u, d) and strange (/ = s) 
quark flavours. In this way one obtains the vacuum masses and scalar densities given in 
Eqs. (0)-©. It is convenient to define the bag constants B f for each flavour separately, as 
the difference between the pressures of physical (m,f — m™ c ) and perturbative (rrif = m^ ) 
vacua, i.e. B f = B(m f) . These quantities are analogous to the bag constants introduced 
in the MIT bag model. For the present choice of parameters, the bag constants are 

B 0u = B od ~ 34.5 MeV/fm 3 , B 0s ~ 208 MeV/fm 3 . (35) 

Note, that for the isosymmetric matter, with equal numbers of u and d quarks, the resulting 
bag constant is B 0u + B od ~ 69 MeV/fm 3 ~ (152 MeV) 4 . 



5 The vacuum will be stable with respect to mass fluctuations if $"(m vac ) < . 
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It is easy to verify that the following thermodynamic identities hold in the model 

p = d li P, p = d 1I P. (36) 
In the mean-field approximation the entropy density 

s = drP = QtPk = 

= -v [nplnnp + (1 - np) In (1 - np) + n p — > Tip (37) 
P 

is the same as for the mixture of ideal gases of quarks and antiquarks. Using the thermo- 
dynamic relations for the free energy density 

F 

T= — = pp + pp-P = e-Ts, (38) 



one can calculate the energy density of the qq matter, e. From Eqs. (B7f), (|3^ ) - (|38|) one has 



e = ^ = e K + ^t + B(m). (39) 

Here the kinetic part, ex, again corresponds to the mixture of quark and antiquark ideal 

gases: 

e K = v E p ( n p + np) ■ (40) 
P 

The same expression for the energy density can be obtained by averaging the energy- 
momentum tensor < T 00 >=< H > /V. According to Eqs. (|29|), (j3~9"D-([40|), the vacuum 
energy density is equal to zero: e vac = B(m vac ) = . 

From the above formulae it is easy to check that the following differential relations 

dP = pdp + pdjj + sdT , (41) 
dT = p dp + ~p dp — s dT , (42) 

hold for any process in a thermally equilibrated system. Using Eqs. (p8[), (f4*2"D one can derive 
the relations 
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p=(pd p + 7dj,-i)r 

= ( P+dp + - l) T + PvdpyJ 7 , (43) 

where p + = (p + p)/2. From this expression it is evident that the pressure vanishes when 
the free energy per particle, F/< AT + iV > = JF/ (p + p), has a minimum in the p — p plane. 
As shall be seen below, at low enough temperatures, a whole line of zero pressure states 
exists in this plane. But only one point on this line, namely the one corresponding to the 
minimum of the free energy per particle should be regarded as a (meta)stable state. At 
T = this state corresponds to the minimum of the energy per particle, 

^ = _f_. (44) 

<N + N> p + p V ; 

Now we discuss the stability of the qq matter. Thermodynamically stable matter may 
exist provided the 2x2 matrix Wdpfip^TW (with pi = p, p 2 = p) is positive 0. If this 
condition is not fulfilled, the system will be unstable due to a spontaneous growth of isother- 
mal density fluctuations (spinodal instability). From Eq. ( fi2]) it is evident that the above 
stability condition holds when the inequalities 

d p p > , (45) 
d p p dpjl > dpp dpjl (46) 

hold simultaneously. As will be shown below, outside the spinodal region, there exists 
a wider region of metastable states where the coexistence of two phases (i = 1, 2) with 
different densities and masses is thermodynamically favourable. This is possible when the 
Gibbs conditions 

pi = p 2 , (47) 
p x = p 2 , (48) 
Pi = P%, (49) 
Ti = T 2 (50) 
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are satisfied. 

In the mixed phase, the condition of the baryon charge conservation can be written as 

X(p 1 -p 1 ) + (l-X)(p 2 -p 2 ) = 3N B /V, (51) 

where A < 1 is the volume fraction of the phase 1 and Nb is the total baryon number of 
the qq system. From Eqs. (ft7|)-(|5"lD one can see that for Nb ^ the pressure of the mixed 
state will in general be A-dependent, i.e. the Maxwellian construction (-P(A) = const at 
T = const) will be violated. The situation is similar to the liquid-gas phase transition 
in nuclear matter with fixed baryonic and electric charges p0j|7|. We would like to stress 
here that by imposing the condition of chemical equilibrium, Eq. (|9|), one normally puts the 
system outside the region of the phase transition. 

In the baryon-free qq matter (< iV — iV >= 0) we have p = p and \i = ~\L (see Eq. (p3|)). 
In this case Eqs. (fJN)— (EH]) become much simpler. In particular, the stability conditions, 
Eqs. pB])-(^o]), are equivalent to the requirement that the isothermal compressibility is 
positive: 

(d p P) T = 2p(d p fi) T >0. (52) 

Here we have used Eqs. ([4lf ), (f45|). 

It is interesting to note that one more instability may occur in a dense baryon-rich 
system due to the presence of strong vector potential. At sufficiently high baryonic density 
a finite qq droplet becomes unstable with respect to the spontaneous creation of qq pairsf]. 
This process becomes possible when the vector field exceeds a certain critical value f4~2fl 

G v \p - p\ > m + m vac . (53) 

When this condition is satisfied, the upper energy levels of the Dirac sea reach the bottom 
of the positive energy continuum. A rough estimate of minimal baryon density necessary to 



This phenomenon is analogous to the spontaneous electron-positron pair production in strong 
electromagnetic fields [41]. 
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produce a qq pair can be obtained by omitting the first term in the r.h.s. of Eq. ( j53|) . In 
the case f — u,d this gives the following estimate for the critical baryon density p% : 

PB = ^(P-P)^f^^6.3p 0) (54) 

where po = 0.17 fm -3 is the equilibrium density of normal nuclear matter. In the above 
estimate we have used our standard choice of the model parameters. At larger Gy the 
critical densities will be lower. 

Before going to the results, let us consider the special case of zero temperature. Taking 
the limit T — > in Eq. fl22|) and using Eq. (p3|) we have 

n p {fJ, R - E p ) = 9(p F - p)6 (p R - m) , (55) 

where 

'a 2 \ 1/3 



^ = v//4-™ 2 = — ( 56 ) 



v 



is the Fermi momentum of quarks. From Eq. fl5"6|) and the analogous expression for anti- 
quarks one gets the explicit formulae for the chemical potentials at T = : 



H = Jm 2 + p 2 F + G v {p-p), (57) 



jj = ^ m 2 + p2 F _ Gv{p _py (58) 
From Eqs. (p0|), fl25D and (|55D we arrive at the following gap equation at zero temperature: 

(59) 



m — m 



Gs 8tt 2 



Two comments are in order here. First, from Eqs. (|57|)-(|5"8D one can see that the condition 
of chemical equilibrium, (|), can not be satisfied at T = . At zero temperature only pure 
baryon (p = 0) or antibaryon (p = 0) matter should be considered as chemically equilibrated. 
Second, according to Eq. (|59|), the quark mass m becomes smaller than m at sufficiently high 
Pf or p F . However, the applicability of the NJL model is questionable at very high densities. 



Indeed, as discussed in Ref. the cutoff momentum A has to exceed all characteristic 



particle momenta. In particular, the inequality A > ma.x(pp,p F ) must not be violated. 
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The explicit formulae for the energy density and the pressure at T = are given by the 
relations: 



e = p,p + pp - P 



8tt 2 




GvPv 



2 



+ B(m) . 



(60) 



The gap equation, Eq. fl59f) , corresponds to the minimum of e at fixed p and p. 



IV. RESULTS 



In this section we present the results of numerical calculations. They were performed 



results are rather sensitive to the particular choice of model parameters. These cases will be 
discussed in more detail. Special attention will be paid to qualitative effects of chemical non- 
equilibrium which, as we shall see later, strongly influence the thermodynamic properties of 
qq matter. 

First, we study the non-strange matter consisting of light (/ = u, d) quarks and an- 
tiquarks. The effects of isotopic asymmetry are disregarded i.e. we assume that p u = pa 
and ~p u = p d . Second, we consider also the ss matter with arbitrary densities of quarks and 
antiquarks. In accordance with the SU(3) flavour symmetry of the interaction Lagrangian, 
the same coupling constants are used for strange and non-strange systems. 

Let us discuss first the properties of the non-strange qq matter at zero temperature. 
As noted above, chemical equilibrium at T = can be realized only at p = (in systems 
with positive baryon charge) . In particular, the baryon-free matter, with equal densities of 
quarks and antiquarks, can not be chemically-equilibrated at T = . By using the formulae 
of Sec. 3 we calculated the quark mass m, chemical potentials p, /Z, the pressure P and 
energy per particle e as functions of densities p and p. The results are given in Figs. 1-7 []. 

7 Unless stated otherwise, for systems with / = u, d we denote by P the total pressure of u and 



by using the parameter set from Eq. @. One should bear in mind that in some cases the 
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Figure 1 shows the density dependence of e, p and m for baryon-free qq matter composed 
of light quarks. The rapid drop of the quark mass reflects the restoration of chiral symmetry 
at high densities. The calculation shows that e(p) has a maximum at a relatively small 
quark density p = Pb ^ Po- At a higher density, p = pa, the energy per particle reaches a 
minimum e = < e(0) = m vac . Therefore, the metastable state of baryon-free qq matter 
is predicted at zero temperature. This state has zero pressure (see below Fig. 3) and its 
existence is possible only in a chemically non-equilibrium system. 

In the considered case the vector density py = and therefore the density dependence 
of e = ef/2pf is determined by the balance between the first kinetic term in Eq. (Q) and 
the bag part B{m) . At p — ► one has m ~ m vac and B(m) ~ ^B"(m vac ) (m — m vac ) 2 . 
The explicit expression for B"(m) follows from Eqs. (p9[), (|3l]). The condition of vacuum 
stability (see the footnote to Eq. fl3~4|)), B"(m v£LC ) > 0, can be easily verified. Substituting 
Pf = Pf = ( 7I " 2 P/) 1 ^ 3 i n t° Eq. (|60|), in the limit ppCm one has 

3 pi B(m) , \ 

10 m 2 Pf y 1 

Minimizing this expression with respect to m, one finds that at low pf the quark mass 
decreases linearly with density: 

2pt 

m ~ m vac - — — . (62) 

B"{m veLC ) 

From Eqs. (|6ll)-(|62|) one can get the approximate formulae 

10 m vac B"(m vac ) 

One can see that at low densities Eqs. (|62])-(|6"3"|) qualitatively reproduce the results shown 
in Fig. 1. Indeed, the second term on the r.h.s. of Eq. (^) is proportional to p 2 ^ 3 . At small p 



d quarks: P = P u + Pd = 2-P u . In the same case p (~p) denotes the total density of light quarks 
(antiquarks), i.e. p = p u + Pd = %Pu ■ When necessary, the single-flavour values are introduced by 
the subscript / . 
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this term dominates and e increases with density reaching a maximum at the point B. From 
Eqs. ([58]) , p3| ) at p = p, T = , one can obtain the relations 

P/2p = pd p e = fi-e, (64) 

which show that P = and /i = e at the extrema of e . 

At p > p( B \ the energy and pressure start to decrease due to the chiral symmetry 
restoration effects. At sufficiently high densities, the quark mass approaches the bare mass. 
In lowest order one can take m ~ m <C Pf and substitute m = m into the expressions for 
P and e given by Eq. (|60"D. This gives the following approximate formulae 

H~p F , e~ -p F + — , (65) 
4 Zp f 

where Bq = B(mo)is the bag constant. From Eqs. d64|)-(|65|) one can see that the pressure 
of a symmetric qq system changes its sign at the point A, where the kinetic (Fermi-motion) 
term and the bag part balance each otherf]: 

Pf A) ^\pfPF-B = 0. (QQ) 

At much higher densities the bag contribution may be neglected and we obtain the limit of 
a massless ideal Fermi-gas: P~|e~|ppiroc p 4 / 3 . 

Substituting B = B 0u from Eq. (|35D we get the estimate pa = 2p| ( A ^ ~ 3.0 po . This is 
close to the density of the metastable state obtained by the numerical calculation, which 
gives the values 

p A ^ 2.93 po, ga = P-a — 0.268 GeV . (67) 

As will be seen below, the value is the absolute minimum of e (p, p) in the whole p — p 
plane. It is interesting to note that the minimum energy value €a is below the vacuum 



"According to these estimates, the origin of the above energy minimum is similar to that suggested 
in the MIT bag model to explain hadron properties in vacuum [ Q . 
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mass, m vac . However, this state is in fact metastable due to the possibility of the annihilation, 
qq — > rnr , n = 2, 3, 

The "chemically-equilibrated" case of baryon-rich matter with p = is considered in 
Fig. 2. A comparison with the preceding figure shows that the quark mass drops with the 
quark density slower than in the baryon-free matter. Although the p-dependence of e is 
rather flat, for the present parameter set the energy per particle has no extrema in the 
considered case. Formally one can use the same procedure to estimate the behaviour of e 
at p — > as was suggested above for the baryon-free matter. The resulting expression may 
be obtained by replacing the last term in Eq. (|63| ) by [Gy — 1/-B"(m vac ) ] P//2. Using this 
expression one can find a fictitious maximum of e at p ~ 6po which is not found in the 
numerical calculation. At such a high density, however, the assumptions m ~ m vac , pp <^ m 
used in the above low density expansion are not valid. Therefore, in the baryon-rich case the 
model predicts a qualitatively different equation of state than that for symmetric qq matter. 

Figure 3 shows the density dependence of the pressure in the u, d matter at T = 0. 
Again, in the baryon-rich case, P(p) is monotonically increasing function. On the other 
hand, the behaviour of the pressure in the symmetric system is non-trivial: there is a 
region of negative pressures at ps < p < Pa- The interval with the negative compressibility 
(d p P < 0) exists between points B and C. One can see from Eqs. (|43|) , (H), that the system 
is mechanically unstable at pb < p < pc- The appearance of regions with negative P and 
d p P are necessary conditions for the existence of bound states and for the first order phase 
transitions, respectively. 

One should bear in mind that the presence or absence of states with negative P and d p P 
in the quark matter with p = is quite sensitive to the choice of the model parameters. 



The similar conclusion was made in Refs. [0,B3|]. By taking small values of the vector 



coupling constant Gy and/or quark bare mass mo, one can find the bound states and phase 
transitions even in baryon-rich systems. This is clear from Fig. 4 which shows the sensitivity 
of the equation of state to the choice of Gy and mo in u, d matter with p = . In the case 
mo = 0, bound states of the baryon-rich matter are possible only if Gy < 0.35 G$- On 
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the other hand, in absence of the vector interaction, i.e. at Gy = 0, the model predicts 
bound states even for bare masses m ~ 5 — 8 MeV, corresponding to correct values of the 
pion decay constant /„. . However, at the realistic choice of the model parameters we do 
not find any phase transitions or bound states in the chemically equilibrated matter. The 
stable quark droplets have been predicted in Ref. |TJ[ within the generalized NJL model. 



As discussed in Ref. pjfl , the prediction is questionable since it was implicitly assumed that 



Gy = 0, m = 0. With the same caution one should regard the results of Refs. fll8l,|20| , l2T 
where the possibility of phase transitions in baryon-rich qq matter was studied within the 
NJL model, either without vector interaction or by assuming unrealistically small values 
of Gy. 

Figures 5(a), (b) show the model predictions for arbitrary qq systems with light quarks 
at T = 0. The contours of equal energy per particle and pressure in the p — p plane are 
given in Fig. 5(a) and 5(b), respectively. Only the vertical and horizontal axes of this 
plane correspond to chemically equilibrated states at zero temperature. One can see that 
symmetric systems with p = p ~ 3 po indeed have the non-trivial minimum of the energy per 
particle and zero pressure. Shaded areas indicate the spinodal regions where the conditions 
(fj5|)-(|4lSD are violated. From Fig. 5(b) one can see that at T = practically all unstable 
states are located inside the contour P = . The presence of a spinodal instability shows 
to the possibility of the phase coexistence in a wider density region. We shall discuss this 
issue in more detail at the end of this section. 

The results for ss systems are presented in Figs. 6, 7. According to Fig. 6, the density 
dependence of m, e and p in the symmetric ss matter is qualitatively the same as for non- 
strange systems. The model predicts a metastable state (point A) with the parameters 

p = p~6.26p , e~ 0.451 GeV. (68) 



Although the assumption m <C pp is not so well-justified as for light quarks, Eq. fl66| ) gives 
the estimate pa — 5.8 po which is close to the numerical result. The noticeably higher 
density of the metastable state in the ss case is due to the higher bag constant of strange 
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quarks. Indeed, from Eq. ( |35"D one can see that B 0s /B 0u ~ 6 . Eq. (|66"D leads to the estimate 
0.5 (B 0s / Bqu) 3 ^ ~ 1.9 for the ratio of densities of the metastable states in strange and non- 
strange systems, rather close to the value obtained from numerical calculations. 

As can be seen from Fig. 7(a), the symmetric ss system with the parameters given 
by Eq. (|68|), has the minimum energy as compared to other points in the p — ~p plane. 
According to Fig. 7(b), the model does not predict states negative P and d p P in baryon- 
rich matter with p = 0. A special study with different vector coupling constants shows that, 
in the latter case, the states with P < exist at Gv/G,s < 0.175. The analogous calculation 
shows that negative values of d p P are possible if Gy/Gs < 0.215. The shaded area again 
shows the region of spinodal instability. For strange qq matter the instability region extends 
to higher densities as compared to the non-strange matter. 

Let us now discuss nonzero temperatures. Eq. (|23|) yields the relations 



P r ( m (RR\ 

-: = I [-, — ), (69) 



V -h&M), (70) 



T 3 V T T 

JL=t (™~Pr 

rp 3 ^ rp ' rp 

connecting the quark (//#) and antiquark (jl R ) reduced chemical potentials with the densities 



p and p . The dimensionless function 



oo , 2 



Jo(x ' y) = i^ l dt — ( 1 2L2 — rrr (71) 

27r z J exp [yx 1 + t 2 — y) + 1 



appears from the integration of the particle occupation numbers in momentum space. 
Using Eqs. (p7|), (p3[), one may write the gap equation in the form 

{d m P K ) T ^ Rr , R = B\m). (72) 

Substituting the expression for Pk from Eq. (p8|) , we have 

B ' (m) -iA^4)+h(ZM- (73) 



mT 2 \T T J \T T 



Here Ii(x,y) differs from I (x,y) by the replacement dt — > dt/^t 2 + x 2 in the r.h.s. 
of Eq. (fnj). The above equation is the analogue of Eq. ( |59"D for the case T / 0. 
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According to Eqs. (§), ([18|)-(p!9D, chemical equilibrium in qq matter is achieved if 
p R = — ~p R f\ Using this condition and Eqs. d69|)-(|70D, (|73|) it is easy to see that the densities 
p and p are not independent in chemically equilibrated systems: the density of antiquarks 
at a given T is fixed by the density of quarks, and vice versa. Figs. 8(a), (b) show the de- 
pendence p (p) at various temperatures for the non-strange and strange matter in chemical 
equilibrium. We see that up to rather high temperatures, T ~ 100 MeV for f = u,d and 
T ~ 150 MeV for / = s, the equilibrium values of p or p are very small. This is especially 
evident for the regions far from the diagonal p = p. 

Figures 9 and 10 represent the mass and pressure isotherms in symmetric qq matter 
with light and strange quarks. These results are obtained by the simultaneously solving 
Eqs. (|69|)-(|70D, (fT3|) with respect to m, p R and JL R . Figs. 9(a), (b) show that the quark 
mass increases with T at a fixed p. Such a behaviour is a specific feature of chemically 
non-equilibrated systems. It can be qualitatively understood, at least in the limit m<T. 
Indeed, from the definition of B(m) it is easy to show that the l.h.s of Eq. (|7^) decreases 
with m. On the other hand, by using Eq. ( p9"| ) one can see that the r.h.s. of the same 
equation increases with p/T 3 . This follows from the fact that Jj(0, y), % = 1, 2 are increasing 
functions of y . This explains the above mentioned behaviour of m (T) . Figs. 9-10 show also 
the density dependence of m and P in the limit of chemical equilibrium, i.e. at p = JI = 0. 
According to Eqs. (|69D , (|73|) , (|27D -(f28D, in this limit m, p and P are fully determined by the 
temperature only. 

As can be seen from Figs. 10(a), (b), at a fixed density the pressure increases with T. 
Similar to the zero temperature case, one can find the spinodal instability regions with nega- 
tive isothermal compressibilities {8 p P)t- The calculations show that the stability condition 



(52]) in symmetric qq matter is violated at temperatures T <T C where 



9 Therefore, the symmetric matter, where pr = p R (see Eqs. (|6^)-([70|) ) , may be chemically 
equilibrated only if p R = Jl R = 0, which in this case is equivalent to p = p = 0. 
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92 MeV, f = s. 



76 MeV, f = u,d, 



(74) 



At lower temperatures the first order phase transition is possible. This means that the 
formation of the mixed phase is thermodynamically favourable. 

The coexisting phases are characterized (for symmetric qq matter) by different densities 
Pi < Pi and quark masses m\ > m 2 , but they have equal chemical potentials (pi = p, 2 = p) 
and pressures. Formally such a situation is possible in the regions of temperature and 
chemical potential where the gap equation has several solutions for the quark mass. Using 
the formulae of Sec. 3, one arrives at the following set of equations for mi, m 2 and p at 



The values of pi (i = 1, 2) are obtained by substituting m = m ; and pr = p into Eq. (^Dj). 
The density points p\ and p 2 are represented by dotted lines in Figs. 10(a), (b). These lines 
constrain the two-phase "binodal" regions shaded at the same plots. At each temperature 
the states with a negative compressibility are situated inside the interval p\ < p < p 2 • 
Similar to the well-known liquid-gas phase transition, the states betwen the spinodal and 
binodal lines are metastable with respect to the heterophase decomposition of qq matter. 

As can be seen from Fig. 10, states with zero pressure exist at sufficiently low temper- 
atures. The maximal temperature when the existence of such states is still possible equals 
T m ~ 49 (67) MeV for non-strange (strange) symmetric matter. At T < T m the states with 
P = are in the metastable density region. Thin lines in Figs. 10(a), (b) represent the den- 
sity dependence of the pressure in chemically equilibrated baryon-free matter. The points 
where this line intersects the pressure isotherms give the corresponding values of P and p 
for the case of chemical equilibrium. These points are located outside of the binodal region. 



fixed T: 



P K (rm,n, T) - B(mi) = P K (m 2 , p, T) - B(m 2 ) , 



< d mi P K (m 1 ,p,T) = B'(m 1 ) 



(75) 



dm 2 P K (m 2 , p,T) = B '(m 2 ) . 
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More detailed information concerning the phase transition in the baryon-free matter is 
given in Fig. 11. The phase coexistence lines in the p — T plane are shown both for non- 
strange and strange baryon-free matter. One can clearly see that states with p — pz — 
do not belong to the phase transition lines. Therefore, the discussed phase trasition has 
essentially non-equilibrium nature. It does not occur in chemically equilibrium systems. 

Figures 12(a), (b) represent the contour lines of e and P for the non-strange qq mixture 
at T = 50 MeV. This temperature is close to the maximum temperature for the existence of 
zero pressure states. The comparison with Fig. 5 shows that the spinodal region occupies a 
smaller part of the p — p plane than in the case T = . From Fig. 8(a) we conclude that at 
such a low temperature the equilibrium quark-ant iquark densitiies lie outside (below) the 
spinodal boundary. An analysis of the binodal boundaries in the p — p plane will be given 
in a future publication. 

V. SUMMARY AND DISCUSSION 

We have demonstrated that the NJL model predicts the existence of metastable states 
of qq matter out of chemical equilibrium. These novel states are most strongly bound for 
symmetric (net baryon-free) systems with equal numbers of quarks and antiquarks. For 
realistic model parameters there exists a region of quark chemical potentials and tempera- 
tures where a first order chiral phase transition is possible. It is shown that for the same 
set of parameters the first order phase transition does not occur in a chemically equilibrated 
system. This conclusion may also be valid for other models of strongly interacting matter. 
Up to now both the chiral and the deconfinement transition in hadronic and quark-gluon 
matter have been studied under the assumption that chemical equilibrium is fulfilled. As 
discussed above, this is most likely not the case in relativistic heavy-ion collisions, at least 
at some stages of matter evolution. 

Let us briefly discuss possible signatures of the phase transitions and of the bound states 
predicted in this paper. Measurements of excitation functions of particle spectra in heavy- 
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ion collisions have been suggested [[L7| to detect a possible critical point in the qq matter 
phase diagram. In principle, the same idea can be applied to the search for these non- 
equilibrium phase transitions. However, an easier task might be to look experimentally for 
the decay of metastable qq droplets in nuclear collisions. It is clear that the formation of 
such droplets would be possible if dense and relatively cold baryon-free matter is created at 
an intermediate stage of the reaction. Such conditions may be realized at high bombarding 
energies, presumably, at RHIC and LHC. Selecting events with unusually high pion multi- 
plicity and low < pt > may increase the probability to observe such multiquark droplets. 
Considerations similar to those suggested in Ref. |^5| show that the subsequent hadroniza- 



tion of droplets may give rise to narrow bumps in hadron rapidity distributions. These 
bumps may be observable in an event-by-event analysis of n, K and spectra. It is also 
obvious, that the yields of BB pairs (B = N, A, S, S, £1 . . .) will be enhanced if they are 
produced through the decay of multiquark-antiquark systems. 

The understanding of the possible decay channels of the predicted metastable states 
is of great importance. In non-strange systems, the strongest decay channel should be 
the annihilation of qq pairs into pions. This channel is energetically open at least at the 
surface of the qq droplet. Indeed, according to Eq. QB7p, the average energy released in the 
annihilation of a qq pair is //+/Z = 2 tA — 3.8 m w . This is sufficient to create up to three pions. 
The largest rate of annihilation is apparently given by the qq — ► 2tt channel. Indeed, the 
one-pion annihilation is suppressed in homogeneous matter due to the energy-momentum 
conservation. On the other hand, the relative probability of the qq — ► 3tt channel should be 
small due to the limited phase space volume available. However, the pion properties in the 
non-equilibrium qq medium have not yet been investigated. It may happen that pion-like 



excitations disappear due to the Mott transition It is also possible that the in-medium 
pion mass is increased as compared to the vacuum value. In both cases the annihilation 
inside the qq droplet will be suppressed. 

The situation in a metastable ss droplet may be even more complicated. It is inter- 
esting that the direct annihilation of the ss pairs into hadrons is suppressed. Indeed, ac- 
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cording to Eq. (|68D , the energy available in annihilation of a ss pair is equal on average 
to /i + /J=2e~ 0.9 GeV. This is lower than the threshold energies of the final hadronic 
states: 0(1020 Me V), KK (990 MeV) and pvr(910MeV). In this case, annihilation is pos- 
sible only for energetic ss pairs from the tails of the thermal distribution. The only open 
channel is ss — > 3 rr . From the analogy to the — > 3 7r decay, which has the partial width 
T^ 3w < 80keV, one can conclude that this channel should be quite weak. 

However, in the deconfined qq matter, the strong decay channels ss — > u u, dd are possi- 
ble, which may reduce significantly the life times of ss droplets. These processes will lead to 
the conversion of ss matter into a system where all three flavours are present. If sufficiently 
long times were available, chemical equilibrium between the flavours would be established. It 
will be interesting to check whether an admixture of qq pairs in strange baryon-rich systems 
is energetically favourable as compared with ordinary stranglets [|ll],|l^]. Detailed calcula- 
tions of the formation probability and life times of metastable qq droplets are possible only 
within a dynamical approach. 

One should investigate also the role of the mesonic degrees of freedom which have been 
disregarded in the present work. It is well known that mesons appear in the NJL model as 
bound states of a quark and an antiquark. The influence of mesonic degrees of freedom in 
chemically equilibrated qq matter has been studied in Ref . |25| . It was shown by using RPA 
that the contribution of mesons may be important at low quark densities. 

In the present calculation the flavour mixing six-fermion interaction was neglected. This 
interaction has been introduced earlier |46|] to reproduce the mass splitting of rj and rj' 



mesons. It has been shown in Ref. p3] that the flavour-mixing interaction only slightly 
changes the EOS of non-strange qq matter at T = 0. It will be interesting to investigate 
the role of this interaction in the case when chemical equilibrium is violated. 
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FIGURE CAPTIONS 

FIG. 1. Energy per particle, chemical potential and constituent quark mass in symmetric 
matter with equal numbers of light quarks and antiquarks at zero temperature, p is 
the total density of u and d quarks. Point A shows the position of the metastable 
state. 

FIG. 2. The same as in Fig. 1, but for asymmetric u, d matter with zero density of anti- 
quarks. 

FIG. 3. Pressure of qq matter with light quarks at zero temperature. The solid (dashed) 
line corresponds to the baryon-free (rich) system. 

FIG. 4. Regions in the m —Gv plane where states with negative pressure (dark shading) and 
negative compressibility (light shading) exist in baryon-rich (p = 0) quark matter 
at T = . Cross corresponds to parameters used in the present model. 

FIG. 5. Contours of equal energy per particle (in GeV, upper part) and equal pressure (in 
GeV/fm 3 , lower part) in qq matter with light quarks at zero temperature. p(p) 
is total density of u and d quarks (antiquarks). Shading indicates the regions of 
spinodal instability. 

FIG. 6. The same as in Fig. 1, but for symmetric ss matter. 

FIG. 7. The same as in Fig. 5, but for ss matter. 

FIG. 8. Quark and antiquark densities, p and p, in chemichally equilibrated qq matter 
(/!= —p) at different temperatures (shown in MeV). Upper and lower part cor- 
respond to the non-strange and strange matter, respectively. 
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FIG. 9. Mass isotherms for non-strange (upper part) and strange (lower part) symmetric 
qq systems. Temperatures are shown in MeV near the corresponding curves. Thin 
solid lines show the density dependence of quark mass in chemically equilibrated 
matter (/i — JI — 0) . 

FIG. 10. Pressure isotherms for baryon-free qq matter. Upper (lower) part represents model 
predictions for non-strange (strange) systems. Thin solid lines show the pressure of 
chemically equilibrated matter. The binodal regions are shaded. 

FIG. 11. Phase coexistence lines in the \i — T plane for symmetric non-strange (solid line) 
and strange (dashed line) qq matter. 

FIG. 12. The same as in Fig. 5, but for the temperature T = 50 MeV. 
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